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Nuclear spin registers in the vicinity of electron spins in solid state systems offer a powerful 
resource to address the challenge of scalability in quantum architectures. We investigate here the 
properties of ^®Si nuclear spins surrounding donor atoms in silicon, and consider the use of such 
spins, combined with the donor nuclear spin, as a quantum register coupled to the donor electron 
spin. We find the coherence of the nearby ^®Si nuclear spins is effectively protected by the presence 
of the donor electron spin, leading to coherence times in the second timescale — over two orders 
of magnitude greater than the coherence times in bulk silicon. We theoretically investigate the use 
of such a register for quantum error correction, including methods to protect nuclear spins from 
the ionisation/neutralisation of the donor, which is necessary for the re-initialisation of the ancillae 
qubits. This provides a route for multi-round quantum error correction using donors in silicon. 


Modular ‘quantum network’ architectures consisting 
of multiple quantum registers connected by interaction 
channels have emerged as a flexible, robust and scalable 
model for quantum computation. Such models typically 
assume high-fidelity operations which can be performed 
locally within the quantum registers (in contrast to po¬ 
tentially lossy channels between them), allowing opera¬ 
tions such as local quantum error correction (QEC) [IHS], 
entanglement purification |4j, and even enhanced quan¬ 
tum sensing mi. This approach is well suited to spins 
of defects in the solid state, such as vacancies in di¬ 
amond [7] or silicon carbide [S], rare-earth dopants in 
various crystals [5] and donors in silicon [TU]. Each of 
these offers a (sparse) environment of nuclear spins, in 
the vicinity of the defect spin, possessing potentially long 
coherence times. This has been explored recently using 
nitrogen-vacancies in diamond, first through the control 
of remote ^®G nuclear spins and later realizing a 

single round of quantum error correction (QEG) [HE]. 

Naturally occurring silicon (''®'*Si) has three stable iso¬ 
topes: ^®Si (92.2 %), ^®Si (4.7 %) and ^®Si (3.1 %), where 
only ^®Si has a non-zero spin (/ = 1/2) and could form 
part of a quantum register. In silicon, much recent focus 
has been on isotopically enriched ^®Si to remove the ^®Si 
spins El: leading to donor electron spin coherence times 
up to 3 seconds m and donor nuclear spin coherence 
times from minutes to hours [Mini- The disadvantage 
of such ^®Si material is that the only additional resource 
for the donor electron spin is the nuclear spin of the donor 
itself. 

Our focus here is on "‘“■‘Si, and in particular the ^®Si 
nuclear spins around the donor. Nuclear spin coherence 
times of ^®Si have been studied in the absence of the 


donor electron (i.e. in bulk NMR m, or using a single 
^®Si atom coupled to a nano-device m) — in such cases 
the nuclear spins can freely flip-flop and the Hahn echo 
T 2 n is limited to around 5 ms. However, the presence of 
the donor electron spin is known to form a ‘frozen core’ 
of nuclear spins around the donor, changing the bath dy¬ 
namics by detuning nuclear spins from their neighbours 
as a result of the spatially varying hyperfine coupling. 
Eor these reasons, one could expect the T 2 n of ^®Si in 
the vicinity of the donor to be significantly longer — an 
indication of this is in the T 2 n of the donor nuclear spin it¬ 
self (strongly detuned from any of the neighbouring ^®Si) 
which was reported to be about 1 second in "‘‘*Si [20] . 

In this Letter, we consider the potential of both the 
donor nuclear spin and local ^®Si spins as a register of 
qubits in silicon, characterising their coherence times and 
examining their use for local QEG. For QEG we consider 
both single-donor approaches (based on single-donor spin 
measurement devices El 1211112) and donor ensemble ap¬ 
proaches (which could form part of hybrid architectures 
with superconducting resonators and qubits |23J [23]). 
In addition to long coherence times, requirements for 
multi-round QEC include qubit manipulation and in par¬ 
ticular the re-initialisation of ancilla qubits. Initialisa¬ 
tion schemes (e.g. by single-spin measurement or optical 
hyperpolarisation [111125]) involve the ionisation of the 
donor, and thus we conclude by examining how to en¬ 
sure a nuclear spin data qubit can be made robust to 
this process. 

We used a float-zone “^-*81 sample doped with 
at a concentration of 6 x 10^® cm“®. Measurements 
were performed at 4.5 K to obtain an electron Ti (> 
5 s) sufficiently long compared to all other experimen- 
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FIG. 1. ENDOR spectra of and nuclear spins 

in "‘‘‘Si at 344.2 mT. (a) Davies ENDOR spectrum of the 
donor in silicon. The peak at 52.475 MHz corresponds 
to a hyperfine of 117.53 MHz, and the linewidth of 60 kHz is 
consistent with previous ENDOR measurements in ""‘Si m- 
RF TT pulse length = 13 fis. (b) Davies ENDOR spectrum of 
^®Si. The hyperfine interaction values with the donor electron 
spin are calculated as twice the shift from the nuclear Zeeman 
frequency of 2.91 MHz, and ranging up to 6 MHz. Inset shows 
a high-resolution spectrum centred around A/2 = 2 MHz, 
showing sub-components of the peaks due to the anisotropy 
of the hyperfine interaction. RF tt pulse length = 50 fis in 
main panel and 1.6 ms in inset. 

tal timescales. Pulsed electron spin resonance (ESR) 
and electron-nuclear double resonance (ENDOR) exper¬ 
iments were realised using a Bruker X-band Elexsys sys¬ 
tem (« 0.3 T, 9.7 GHz). The magnetic field was set 
parallel to the [001] Si crystal axis, where the electron 
spin coherence time T 2 e is maximized (« 0.5 ms [26)1. 
Dynamical decoupling sequences applied to the nuclear 
spins were synthesised directly from an arbitrary wave¬ 
form generator (Agilent 81180). 

We begin by characterising the and ^®Si nuclear 
spins through Davies ENDOR spectroscopy [mUH]. The 
^‘P donor nuclear spin has a well-known gyromagnetic 
ratio of 17.23 MHz/T and a hyperfine interaction with 
the donor electron spin of 117.53 MHz [29]. ^®Si spins 
in the bath around the donor have a gyromagnetic ra¬ 
tio of 8.46 MHz/T and hyperfine couplings to the donor 
electron spin of up to 6 MHz. An in-depth study of all 
couplings and related sites can be found in EOj. Spec¬ 
tral overlapping makes weakly coupled ^®Si more difficult 
to distinguish experimentally — for these, the hyperfine 
interactions can be simulated using the Kohn-Luttinger 
model of the electron wavefunction (see Supplementary 
Material). 

We then measure the coherence times ( 7211 ) for these 
various nuclear spins, based on the approach of coherent 
state transfer from the donor electron spins to the nu¬ 
clear spin, and back again at some later time, as shown 
in Fig. [^a) [101. Microwave pulses on the ESR transi¬ 
tions must be selective on a particular nuclear spin state, 
and thus have a bandwidth significantly less than the 
relevant hyperfine coupling. This is trivial in the case 



FIG. 2. Nuclear spin coherence times of ^‘P and ®®Si. 

(a) Left: energy diagram for the donor electron spin coupled 
to a spin-1/2 nuclear spin. Right: nuclear spin coherence mea¬ 
surement sequence taken from Ref pro). The blocks defined by 
the dashed brackets move together when r is varied, (b) 
nuclear spin coherence decay. The signal shown is the mag¬ 
nitude of the ESR in-phase and quadrature detection, hence 
the fit (red) is constrained to decay to zero as the noise is 
always positive, (c) Coherence decay for a ®®Si nuclear spin 
with A = 4.03 MHz. 


of the donor nuclear spin, however for spins this re¬ 
quired microwave pulse lengths of 0.5 /js. Figure [^b) 
shows the nuclear spin coherence decay is observed for 
with a resulting decay time T 2 n = 1-1 ± 0.1 s. A 
comparable coherence time of 1.22 ± 0.03 s was mea¬ 
sured for a ^®Si nuclear spin with hyperfine coupling 
A = 4.03 MHz (Fig. [^c)), notably over 200 times longer 
than in bulk natural silicon. In both cases, the decay 
followed a stretched exponential function exp (—'r/T 2 n)" 
with stretch factor n around 2, typical of decoherence 
from spectral diffusion in ""‘Si |5T] . 

The ^®Si nuclear spin coherence time was found to de¬ 
pend strongly on the hyperfine coupling to the donor 
electron spin, as shown in Fig. [^a). For the strongest 
hyperfine coupling (A > 3 MHz), the coherence time sat¬ 
urates at « 1.3 s, and then decreases with weaker cou¬ 
pling, towards the bulk NMR value of 5 ms [TH]. Two 
decoherence mechanisms can be considered in this case. 
The first is due to flip-flop of ^®Si spin pairs very far 
from the donor, and thus far from the measured ^®Si 
spin. The distant ^®Si spin pairs have negligible hyper¬ 
fine interaction with the donor electron spin (compared 
to the dipolar coupling within the pair), and are therefore 
not detuned from one another, allowing for flip-flops [32] . 
Their small coupling with the measured nuclear spin is 
compensated by the very large number of pairs involved 
in the process (« 10® [33]). A second process, recently 
proposed |33j . involves only a few spin pairs that are 
much closer to the donor and are located at lattice sites 
that are equivalent by symmetry. Such pairs have equal 
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FIG. 3. coherence time as a function of hyperfine 
coupling and dynamical decoupling (DD) (a) The 2®Si 
coherence times, T 2 n, vary with the strength of the hyperfine 
coupling to the donor electron (and thns, indirectly, as a fnnc- 
tion of the distance between the two). The line in black is only 
a guide to the eye, limited at low coupling to the bulk NMR 
value (5 ms) and at large coupling to « 1.3 s. (b) Measured 
T 2 u times under various DD sequences for a specific ^®Si site 
with coupling A = 2.23 MHz. CPMG, tt-WAHUHA, and XY- 
4 all offer identical protection of the nuclear spin coherence 
(in proportion to the number of refocusing pulses), showing 
that indirect flip-flops in the environment of the measured 
^®Si spin are responsible for decoherence. 

coupling to the electron spin, hence there is no detuning 
within the pair. Both of these mechanisms induce deco¬ 
herence to the measured nuclear spin via the Ising (ZZ) 
interaction, also termed an indirect flip-flop process as it 
is due to flip-flops of neighbours [31]. We go on to use 
dynamical decoupling (DD) to further investigate these 
decoherence mechanisms. 

Dynamical decoupling has been used extensively in dif¬ 
ferent contexts ranging from extending coherence times 
OUT] and performing spectroscopy [HIES] to probing 
quantum interactions [36l|37]. In Fig. [^b), ^®Si nuclear 
spins at a specific site {A = 2.23 MHz) are subject to 
different types of DD sequence: CPMG [Ref [3B]], XY- 
4 [Ref jUj] and a modified version of WAHUHA [Ref [lO]]. 
CPMG consists of a train of tt pulses that refocuses 
ZZ interactions between spins. Our experiments show 
that under CPMG r 2 n increases linearly with the num¬ 
ber of TT pulses, up to 3.7 ± 0.2 s (for eight tt pulses). 
This improvement provides an additional evidence that 
indirect flip-flops are the likely source of decoherence. 
This can be further tested using the WAHUHA sequence: 

mod¬ 
ified here (called tt- WAHUHA) to also include tt pulses to 
allow for refocusing of inhomogeneous broadening {T 2 ) ■ 
By alternating the rotation axis of the | pulses, this re¬ 
focuses the dipolar interaction between the measured nu¬ 
clear spin and any equivalent pair. By comparison with 
the results from CPMG, it can be seen that this sequence 
does not improve the nuclear spin coherence beyond what 
would be expected from its five tt pulses, which eliminates 


the possibility of a decoherence mechanism due to direct 
flip-flops. Finally, XY-4, which has four tt pulses with al¬ 
ternating rotation axes, is applied to check for any effect 
from pulse errors, and unlike CPMG is a universal DD 
protocol required for use in general qubit applications. 
In summary, the coherence of both ^^P and ^®Si nuclear 
spins have been measured to be in the order of seconds, 
and can be extended using DD sequences. 

The long coherence times measured above demonstrate 
that nuclear spins near the donor could be used as a 
quantum register, however, applications such as quan¬ 
tum error correction require the ability to repetitively 
initialise the states of ancilla qubits. Even at low tem¬ 
peratures (< 100 mK) and high magnetic fields (> 1 T), 
the nuclear spins are in a fairly mixed state in thermal 
equilibrium, however, the polarisation of the donor elec¬ 
tron spin can be transferred to the nuclear spins, fol¬ 
lowing the same methods used in the ENDOR experi¬ 
ments above. Two methods to polarise donor electron 
spins quickly and on-demand include i) the use of spin- 
selective donor ionisation through the use of the bound- 
exciton IR transition (applicable in both ensembles and 
single spins) US] ST] |33| ; and ii) the measurement of a 
single donor spin coupled to a single electron transistor 
(SET) |43|. In the first case, laser excitation (at around 
1078 nm for ^^P) causes only donors of a defined spin 
orientation to be ionised, which is followed by a subse¬ 
quent capture of an electron in a random spin state. This 
can achieve full donor electron spin polarisation on the 
tens of millisecond timescale (depending on laser power). 
Although the strain caused by the isotopic variation in 
natural silicon leads to a broadening of the donor-bound 
exciton linewidth, the electron spin can still be resolved 
at modest fields (see Fig. |^a)). In the second case, the 
timescales are set by tunneling rates between the donor 
and the SET, which give a measure/reset time of order 
1 ms. 

Both of these spin initialisation methods rely on ioni¬ 
sation of the donor, which impacts the coherence of any 
coupled nuclear spins in two distinct ways. First, while 
the donor is ionised there is no longer a ‘frozen core’ of 
protected nuclear spins and so the flip-flops in the nuclear 
spin bath limit r 2 n to the 5 ms timescale [T^. During 
such periods, DD sequences similar to WAHUHA can 
be applied to suppress the dipolar interaction between 
the spins, as was already demonstrated using NMR in 
Ref. [44] where the ^®Si nuclear spin coherence was ex¬ 
tended up to 20 s. A second issue arises from the in¬ 
herent uncertainty in the precise timing of the ionisa¬ 
tion/neutralisation of the donor, as this imparts a ran¬ 
dom phase on the nuclear spin related to the strength of 
its hyperfine coupling to the donor electron. If the nu¬ 
clear spin state is an eigenstate, it is rather insensitive to 
the donor ionization, as evidenced by both optical and 
electrical ionisation experiments naEH, however while 
it is in a superposition state one can expect the ran- 
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dom timings of the donor electron removal/re-capture to 
lead to decoherence. Notably, this decoherence process 
is also observed in nuclear spins near NV centres in di¬ 
amond where prolonged measurement of the NV centre 
can cause it to randomly change its charge state [7]. 

One solution is to use nuclear spins whose coupling 
to the donor electron spin is much weaker than the in¬ 
verse of the ionisation time uncertainty, but this would 
require using ^®Si with hyperfine values <?; 1 kHz, which 
in turn have short coherence times and whose condi¬ 
tional operations through the donor electron spin would 
be slow. We hence suggest protecting the nuclear spin 
coherence by applying DD on the electron spin, at times 
when ionisation/neutralisation of the donor is expected. 
The hyperfine interaction can thus be effectively turned 
off on-demand, assuming that the pulses are applied at 
a repetition rate much faster than the hyperfine cou¬ 
pling strength (see Supplementary Material for deriva¬ 
tion). Critically, the hyperpolarisation control (in the 
form of laser or voltage pulses) must be synchronised 
with the DD pulses in order to work effectively, as il¬ 
lustrated in Fig. I^b). Following this protocol, the elec¬ 
tron spin state can be reinitialised while the coherence of 
(weakly coupled) nuclear spins remains unperturbed (see 
Fig. Qc)). Finally, this DD method could have further 
applications, such as protecting the nuclear spins from 
Tie relaxation of the electron spin (similar to Ref |45) 1. 

Further considerations (see Supplementary Material) 
for the implementation of a quantum register based on 
^®Si weakly coupled to the donor include i) the effect 
of anisotropy in the hyperfine coupling, and ii) shot-to- 
shot fluctuations in the state of the nuclear spin bath 
(manifest as a ESR linewidth of « 8 MHz Ref [2T] ). The 
former could lead to undesired nuclear spin flips as a 
result of DD applied to the donor electron spin, and can 
be mitigated by increasing the magnetic field strength. 
The latter shifts the ESR frequency, however, conditional 
operations on the electron spin controlled by the nuclear 
spin can still be performed through the use of frequency 
comb methods (similar to Mims ENDOR [M]). 

In conclusion, we have considered the suitability of ^®Si 
nuclear spins around a donor electron spin as a quantum 
register, and measured their coherence times to be in 
the seconds timescale and a function of their hyperhne 
coupling to the donor. These could be harnessed to per¬ 
form, for example, a three-qubit QEC protocol using the 
donor nuclear spin and one strongly coupled ^®Si as an- 
cillae, and one weakly-coupled ^®Si for the data qubit. 
Combined with recent measurements which show that 
bismuth donor electron spin coherence times can reach 
a second in natural silicon 07], these results indicate 
that isotopically enriched ^®Si may not be a panacea for 
silicon-based qubits, and the more abundant and easily 
accessible variant may bring benefits for some applica¬ 
tions. Although more technically complex, there may 
also be merits in incorporating ^®Si in the vicinity of the 
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FIG. 4. Resetting the donor electron spins whilst pre¬ 
serving nuclear coherence (a) Donor bound exciton (D°X) 
energy diagram and measurement of electron spin hyperpo¬ 
larisation P = polarized ^ ^onors in "“-‘Si, 

where S'Echo polarized and 5*Echo thermal are the electron spin 
echo signal intensities with and without illumination. ESR 
frequency, / = 9.7 GHz, 349 mT, T = 4.5 K. The actual 
spin polarisation might be somewhat smaller as the enhance¬ 
ment observed could include a contribution for donor nuclear 
spin polarisation, due to cross-relaxation, (b) Sequence for 
protecting a weakly coupled ^®Si nuclear spin coherence dur¬ 
ing donor spin hyperpolarisation by spin-dependent tunneling 
(voltage pulses, low level = “read” and high level = “load” 
according to Ref [mils]) or spin-selective optical ionisation 
(laser pulses). DD on the donor electron spin (ESR) is syn¬ 
chronised with the laser/voltage pulses in order not to disturb 
the electron spin polarisation process. DD on the nuclear spin 
(NMR) is a WAHUHA-like sequence with 7r/2 pulses to refo¬ 
cus the dipolar interaction, protecting the nuclear spin coher¬ 
ence from flip-flops when the donor is ionised, (c) Simulation 
of the sequence in (b) in the case of spin-dependent tunneling, 
showing the evolution of the donor electron spin and charge 
states, and the ^®Si nuclear spin coherence with and without 
DD on the electron spin. Electron spin populations are plot¬ 
ted after every other ESR tt pulse. Simulation parameters: 
donor ionisation and capture characteristic times are 295 fis 
and 33 /is, respectively (taken from Ref |21p. rr-pulse de¬ 
coupling rate is 5 MHz for a hyperfine interaction strength of 
0.1 MHz. Spin-selective optical ionisation (laser pulses) shows 
similar evolution but on longer timescales (10-100 ms). 


donor (e.g. through co-implantation), in material which 
is otherwise isotopically enriched. 
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